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Abstract 

In this paper we derive a necessary optimality condition for a local opti- 
mal solution of some control problems. These optimal control problems are 
governed by a semi-linear Vettsel boundary value problem of a linear elliptic 
equation. The control is applied to the state equation via the boundary and 
a functional of the control together with the solution of the state equation 
under such a control will be minimized. A constrain on the solution of the 
state equation is also considered. 
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1 Introduction 

In this paper we discuss the necessary optimality conditions for a class of optimal 
boundary control problems governed by a linear elliptic partial differential equation 
with nonlinear Vettsel boundary condition. We formulate the problem first. 

Let J : C a (dfl) — > M. be the objective functional defined by 

J(u)= / f(x,y u (x)) dx+ / g(s,y u (s),u(s)) ds (1) 
Jn Jan 
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where / : Q x R — ► R and g:9flxIxR^I are of class C 1 and y u = G{u) is the 
solution of the state equation 

-Ay + y = in SI, B(x,y,u) = on (2) 

corresponding to w. The boundary condition B(x,y,u) = in the state equation 
will be defined in details in Section 2. Let Z7 be a bounded set of functions, called 
the set allowable controls, in C a (dfl) and F : C 2,a (Q) ^Ibea constrain functional 
on the state y u given by 

F(u) = / a(x,y u (x)) dx + / b{s,y u (s)) ds (3) 
7n Jan 

where a : f2 x R — > R and 6 : <9fi x R — > R are both of class C 1 . The function spaces 
C a (dQ) and C 2 ' a (fi) will also be introduced in details in Section 2. The control 
problem is formulated as follows: 

Minimize J(u) 

u(x) E U (4) 
F{y u ) = 

When the state equation is an elliptic equation with a traditional boundary con- 
dition, i.e. either a Dirichlet or a Neumann boundary condition, or a combination 
of the two in the form of a general oblique boundary condition, this problem has 
been well studied. A survey of those results is given in [3]. For first and second 
order necessary optimality conditions for state equations with Neumann boundary 
conditions we refer to [I] and [6]. 

However problem (j3j) has not been studied so far when the boundary condition 
is of Vettsel type. A Vettsel boundary condition consists of not only the unknown 
function and its first order derivatives but also the second order tangential derivatives 
of the unknown function. It has been shown in [TT] that Venttsel boundary condition 
is the most general admissible boundary condition for second order elliptic operators. 

Such a boundary condition can also be found in many engineering problems and we 
refer [TJ and [8] for details. A simple example of it is the problem of heat conduction 
in a medium enclosed by a thin skin and the conductivities of the medium and the 
surrounding skin are significantly different, see [3|. Generally speaking, all physical 
phenomena involving a diffusion process along the boundary will give rise to a Vettsel 
type boundary condition. 

In the following, we will first discuss the existence theory of the state equation 
in Section 2. In Section 3 We will establish the differentiability of all functionals 
associated to our problem and find the derivatives of them. In Section 4, we will 
state and prove the main theorem of this paper. Finally in Section 5 we will make 
some comments on work in progress and further development. 
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2 State Equation 



Let Q be a bounded open subset of R n with a C 3 boundary dQ. Let : dQ x R x 
R —>■ R be a C 1 function. Given a function w 6 C Q (<9f2) we consider the following 
boundary value problem: 

- Ay + y = in fi, 

A m y + d„y = (p(-,y,u) on <9fi (5) 

The boundary condition in (jHJ) is a special case of the general Venttsel boundary 
condition. The definition of such a boundary condition is given as follows. Let 
v = (u 1 , . . . , v n ) be the inward unit normal vector field on dVt. Then the inward 
normal derivative of u, denoted by d u y, is defined by 

d v y = Dy-v 

where Dy is the gradient vector. Now we define the tangential differential operators. 
Let {c lk } nxn be the matrix whose entries are given by 

c ik = 6 ik _ v i y k 

where 5 th is the Kronecker symbol. Then the first and the second order tangential 
differential operators are then defined by 

di = c lk D k , dij = didj, i, j, = 1, • ■ ■ , n, 

hence the tangential gradient operator is defined by 

d = (<9i, • • • ,d n ). 

In particular the Laplace-Beltrami operator on the boundary manifold is then de- 
fined by 

A 9n = didi. 

Note that the second order tangential derivatives so defined are not symmetric in 
general. 

In order to understand the optimal control problem, the fundamental issues is the 
existence and uniqueness of solutions of the state problem ([5]), as well as the contin- 
uous dependence of the solutions upon the input u. The existence and uniqueness of 
solutions of linear problems have been studied in [8] . Similar results for quasi-linear 
equations with quasi-linear boundary conditions are covered in [9]. Our problem ([5]) 
is a linear equation with a semi-linear boundary condition, so the solvability and 
uniqueness can be deduced from the general frame work of [9] . In order to make the 
article more readable, without having to verify those complicated general structure 
conditions stated in [9] , we prove the existence and uniqueness using only the results 
of 0. 
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We seek classical solutions in the Holder space C 2 ' a (Q). For a non- negative integer 
k and a number < a < 1 the general Holder space C k,a (Q) is the Banach space 
whose norm is defined by 

\y\k,a-,si = \y\k,o-,ci + [D k y] a -,n 
where D k y denotes the fcth order partial derivatives of y, 

k 

\y\k,a-,n = supV" \D % y(x)\ 

and 

[nfe i \D k y(xi) -D k y(x 2 )\ 
[D y\ a] n = sup . . 

Notice that every C k ' a (dQ) function can always be extended to a C k ' a (Q) function 
and, on the other hand, every C k ' a (Cl) function can be restricted on the boundary 
to produce a C k ' a (dQ) function. Notice also that both conversions can be carried 
out in a manner that preserve the norm, i.e. the corresponding C k,a (dfl) norm and 
C k,a (Cl) are equivalent. Based on such an observation we will not distinguish the 
spaces C k ' a (dn) and C k > a (Cl). 

For the existence and uniqueness we make the following assumptions. The definition 
of the allowable set of control implies that there is a constant Mi such that 



\u\ ,a-,n < Mi. 



We assume that ip satisfies 



1. For x G fi, y G K. and |w|o ; n < Mi, there is a positive constant cq such that 

dip 

Q-(x,y,u) > c > 0. (6) 

2. For i6fl, \y\o-n < Mo and \u\o-n < Mi, there is a constant M 2 depending on 
Mo and Mi such that 

\<p(x,y,u)\ , \D(p(x,y,u)\ , \D 2 <p(x,y,u)\, \D 3 (p(x,y,u)\ < M 2 . (7) 

Let us now consider the linear problem 

- Ay + y = mfl, A dn y + d u y = ay + h on dQ (8) 
where a, h G C a (dQ) are given functions such that 

a(x) > c > 0. (9) 
We may assume that the c in ([9]) and ([6]) are the same. 

As a starting point we quote Lemma 1.1, Theorem 1.5 and Theorem 1.6 of [8] here. 
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Lemma 2.1 If Xq G dQ is a maximum point of y G C 2 (Q) then at x we have 
AanV < and d u y < so that 

&dny + d v y < 0. (10) 



Theorem 2.2 Suppose that Q is a C 2,a domain and (TJ|) holds. If y G C 2,a is a 
solution of (EJ) then 

\y\2,a;Q < C(Mo ; fi + Mo.ojfln) (U) 

where C only depends on the geometry ofdQ, cq, |a|o, Q ;n andn. 

Theorem 2.3 If Q is a C 2,a domain then the boundary value problem (E[) has a 
unique C 2,a solution for every h G C 2,a (dQ). 

To handle nonlinear problems we frequently need the following fact: 

Lemma 2.4 Suppose ip : Q x W k ^ M. is a C 1 function satisfying 

M,|£ty|<M 2 

for a constant M 2 . Ifui(x) G C a (Cl) for i = 1, 2, . . . , k then ip(x, Ui(x), 
C a (tt) and 

< M 2 ^1 + d l - a + £ | 
where d is the diameter offl. 



■ ■ -,u k (x)) G 
(12) 



Proof. Obviously 

Also by the assumption on if> we have 

(xi,Mi(xi), . . . ,U k (xi)) - 1p(x 2 ,U 1 (x 2 ), ■ ■ .,Uk{x 2 ))\ 



sup 



< sup 

X\ ,X2d£l 



dip 



dx 



\Xi - x 2 



\xi — x 2 \ a 

k 



8=1 



dip 



dui 



\Uj(xi) -Uj(x 2 ) \ 
\x\ — x 2 \ a 



< M 2 [d l - a + Yy 



i=l 



Adding these up yields (TT21) . 

Now we turn to our nonlinear problem (JSJ). Suppose that for each u 6 (/ the problem 
(loT) has a classical solution y. We show first that the C°(Q) norm of y is bounded 
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by the C°(Q) norm of u and therefore it is bounded as U is a bounded set in C°(Q). 
Let M = \y\o-,ci = sup n \y\ > 0. Without loss of generality we assume M = y(x ) 
for some Xq G 0. By the weaA: maximum principle we have £0 G <9f2. Then Lemma 
12.11 implies 

> A a ny(x ) + d u y(x ) = ip{x , M , u(x )) 

= M —{x ,tM ,u{x ))dt + if{xo,0,u{x )) 
Jo dy 

> c M + ip(x o ,0,u(x Q )). 
Thus by (JSJ) and ([7j) we have the desired bound for y: 

M <-sup|^(-,0,«(-))|<-M 2 . (13) 
Co n c 

Next we show that the C 2,a (Q) norm of y is also bounded. For this purpose we put 
h(x) = <p(x,y(x),u(x)). For convenience we may assume, from now on, that d > 1 
so d}~ a < d. Then, by Lemma [2.41 we have 

|/i|o,«i«i < M 2 (l + d + M l + \y\ (14) 
By inserting (1131) and (j!4p into (II ip we obtain 

\y\2 >a; n < C(M + M 2 (l + d + Mi + Mo,a;n)) = Ci + C 2 |y|o,a;n (15) 

for some constants C\ and C 2 . Now we recall the well known interpolation inequality 
(Lemma 6.35 of [7j): 

|z/|o,a ; n < C|y|o : n + e|2/|2,a;fi < —C 3 M 2 + e\y\ 2 , a -.n (16) 

Co 

for some constant C = C(e,Q). By choosing e = 1/(2C 2 ) we have proved the 
important a priori estimate for the solutions of (151): 



Theorem 2.5 Assume that ip satisfies (EP and (0), and u £ U. If y is a solution 
of (EP i/ien we /jat> e 

|2/|2,a;n < C 

for a constant C depending on Mi, M 2; Q and n. 



With u fixed, for all z G C 2,a (dQ) we define the operator T by letting y = Tz be 
the unique solution in C 2 ' a (dQ) of the linear problem 

- Ay + y = in Q, 

Adtiy + d u y = ip{-,z,u) on d£l. (17) 

The unique solvability of this linear problem is guaranteed by Theorem 12.31 By 
applying the Leray-Schauder fixed point theorem together with Theorem 12.51 to the 
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operator T we obtain the existence of a fixed point of T which is then a solution of 
our state equation fl5J). 

To see whether the solution of (jHJ) is unique we suppose that 2/1 and 2/2 are two 
solutions. Then the difference y = 2/1 — y2 will satisfy 

- Ay + y = in O, 

^dQV + d v y = (p{-,yi,u) - (p(-,y 2 ,u) on 90. (18) 

If 7^ 2/2 we may assume M = 2/1(^0) — 2/2(^0) — su Pn — 2/2 1 > 0. The weak 
maximum principle implies that xq G dQ. Then by Lemma [2.11 at xq we have 

> A 9Q y + = v?(-, 2/1, u) - ¥>(•, 2/2, w) > 0, 
a contradiction. In summary we have 

Theorem 2.6 Assume that (p satisfies (TJ|) and |7p, and u £ U. Then the state 
equation (Tj5|) /ias a unique solution y G C 2,a (f2). 



3 Differentiability 

Given u G C a (<9f2) we denote by y u the solution of the boundary value problem (jSJ). 
This correspondence defines a mapping G : C a (dQ) — > C 2 ' Q (f2) by t/ u = £?(«). 

Theorem 3.1 The mapping y = G(u) is Frechet differentiable. Let G'(u) be the 
Frechet derivative of G and z = (G'(u), v) where v G C a (dQ). Then z is the unique 
solution of the boundary value problem 

- Az + z = in n, (19) 

d d 

A dn z + d v z = — {-,y,u)z + —{-,y,u)v on dtt. (20) 
ay ou 

Proof. We first prove that G is Gateaux-differentiable and calculate the G- 
derivative dG{u). Let v G C a (dQ) and consider y t = G(u + tv) and y = G(u). 
It follows that 

- Ay t + y t = in Q, 

A dn y t + d v y t = <p(-,y u u + tv) on dtt (21) 

and 

- Ay + y = in Q, 

A 9a y + d u y = ip(-,y,u) on <91X (22) 
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By subtracting fl22|) from (1211 we see that w t = yt — y satisfies 

— Aw t + w t = in Q, 

A d ciWt + d„w t = ip{-,yt,u + tv)-ip(-,y,u) on <9fi. (23) 

We can assume that t is bounded, say \t\ < 1. Theorem 12.51 guarantees that the 
C 2 ' a (Q) norms of both y and yt are uniformly bounded. Once u and v are chosen 
their norms are also independent of t. We now want to show that the norm |wt|2, a; n 
is also uniformly bounded with respect to \t\ < 1. For this purpose we consider the 
composite function h(x) = ip(x,y t (x),u(x) + tv(x)) — <p(x, y(x), u(x)). By Lemma 
12.41 again we have 

\h\ ,a;a := \h\ -,n + [h)a-,n < C (24) 
for a constant C that is independent of t. Using Theorem 12.21 we fist obtain 

\w t \ 2 ,a;n < Ci(K|o ; n + \h\ ,a;dn) < Ci(ki|o ; n + C) (25) 

Then the problem is reduced to the estimation of |wt|o ; n- The previous argument 
does not work here because equation (123]) does not have the w t term. However this 
situation is covered by Lemma 1.4 of [8] which gives 

Klo ; n < C 2 sup \<f(x,y t (x),u(x) + tv(x)) - ip(x, y(x), u(x))\ < 2C 2 M 2 (26) 

where the constant C 2 is independent of t. A substitution of ( l26l) into ( |25l) produces 

\w t \ 2 ,a;n < C 3 (27) 

for a constant C3 independent of t, that is, |iUt|2,a;n is uniformly bounded. Therefore, 
up to a subsequence, the following limits exist in C 2,a (Q): 

limu> t = w and lim% = \im(w t + y) = w + y. 

t->o t-*o t-»o v 



By taking limit in (1231) as t — > we have 

— Aw + w = in Q, 

A dn w + d u w = <p(-,y + w,u) - ip(-,y,u) on dQ. (28) 

Since ip is increasing in y variable, the only solution satisfying (128]) is w = 0. From 
this we conclude lim^o Vt — V- 

Now consider z t = w t /t. Dividing ( 123]) by t yields 

— Az t + z t = in Q, 

A dn z t + d u z t = a t z t + (3 t v on dVt (29) 

where 



l -^—(x,Ty t + (l — T)y,u + tv)dT and f3t = -—(x,y,u + Ttv)dr. 

Jo oy Jo du 



Obviously a t G C 2,a (Q). It follows from Lemma [2.41 that 

|a|o,a;n < M 2 (d + \y\o, a; n + \yt\o,a-,n + \u\o, a -n + \v\o, a; n) < C 4 

for a constant C4. Notice also that at > cq and \(3t\ < M2 and hence Theorem 12.21 
implies 

N2,a;Q < CE)(No;n + |A^|o,a ; n) < C 5 (|^|o;fi + C 6 ) 

for some constants C5 and Cg. Finally the uniform bound for |zt|2,a;fi comes from 
the estimate 

c 

which is a consequence of Lemma 12.11 In summary we have 

Iztkojn < ^7 (30) 

for a constant C7 that is independent of t. This implies that, up to a subsequence, 
z t converges to a function z in C 2 ' a (f2) as £ — > and 

lim a t = ^(x,y,u) and lim f3 t = ^-(x,y,u). 
t^o ay t-*o du 

Taking limit in (129]) gives 



- Az + z = in Q, 

d d 

A dn z + d u z = —(-,y,u)z + —(-,y,u)v on dtt (31) 
ay ou 

which means that z = (dG(u),v) is the solution of f)19p and (T2U1) . 

The uniqueness of 21 is guaranteed by Theorem 12.31 as f[3"Tl) is a linear equation. 

Next we examine the continuity of dG. Notice that dG{u) G C[C a {dVt) ) C 2,a {pl)) 
and 

= sup \(dG(u),v)\ 2 ,a;n- 

\\v\\=l 

Therefore to prove the continuity of dG(u) is to prove that as u% — > u in C a (Cl) 
\\dG{u\) — dG{u)\\ = sup \(dG(ui),v) — (dG(u),v)\2 t a-,n — ► 0. 

For any v G C a (dVL) with ||t>|| = \v\ a] a = 1 consider z\ = (dG(ui),v) and z = 
(dG(u),v). Then we know that W\ — Z\ — z is a solution of 

— Awx + w\ = in f2, 

a o cfy? dep. , 

Aen^i + d^wx = —(•,?/, «i)zi + — (■,?/, 
ay aw 

d d 
9 



All we need to show is that W\ — > in C 2,a (Q) uniformly with respect to \v\ a .Q = 1, 
as U\ — > u in C a (0). To this end we rewrite equation (13"21) in the form 

— Awi + wi = in Q, 

AdnWi + d v W\ = crwi + 7 on dQ (33) 

where 

° = d~( x,y,Ul ) 

and 

I ' ®V < \ dip / A /cV/ \ dip . A 

7 = ^(-,y,«i) - ^ (•,*,«) J*+(^0.y.«O - «) J «■ 

If we put 

A= (x,y,u + t(u\ — u)) dt and B= I ——(x,y,u + t(ui — u))dt 

Jo dudy Jo du 2 

then 7 can be written as 

7 = (Az + - ux). 

From the assumption on ip we know that (Az + Bv) G C a (n) and hence 

By Theorem 12.51 we then have 

ki| 2 ,a;n < Cgd^ilojn + ItUo) < Cgdwilojn + C 8 |wi| a; n) < C 9 (l + C 8 )|iOi| a;n . 

Now the continuity of dG{u) follows because |u>i| a; n — > uniformly with respect to 
\v\ a -,u = 1- Finally, since G(u) is continuously Gateaux differentiable, we conclude 
that G(u) is also Frechet differentiable and that the Frechet derivative G'(u) is equal 
to dG(u). 

Now we are in the position to establish the differentiability of the objective functional 
J{u). 

Theorem 3.2 The functional J is Frechet differentiable and for every u, v G C a (dQ) 
and y = G{u) we have 



(J'(u),v) 
where w is the solution of 



Ml 



dg , \ dip , ' 
-(s,y,u)--(s,y,u)w 



v ds 



df 

-Aw + w=-^-(-,y) in tt, (34) 

dy 

A dn w + d u w = ^-(-,y,u)w-^-(-,y,u) on dtt. (35) 
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Proof. Define 

H(y,u)= / f(x,y(x)) dx + / g(s,y(s),u(s)) ds. 
Jn Jan 

It follows that 

J{u) = H(G(u),u). 
It is obvious that H is differentiable and for every y and u we have 

(^(y,u),y) = [ ^f(x,y(x))ydx + [ ^-(s,y(s),u(s))y ds 
oy Jn oy J an dy 

and 

dH f dg 

(-^-(y,u),u) = — (s,y(s),u(s))uds. 
du J dn du 

By the chain rule we have 

PjTT f)M 

(J'(u),v) = (—(y,u)G\u) + -^(y,u),v) 

which then gives 

(J'(u),v) = f tt(x, y(x))G'{u)v dx + f ^-(s,y(s),u(s))G'(u)v ds 
Jn oy J d n oy 



Jan ou 



df f dg 

— (x, y(x))z(x) dx + / — (s, y(s), u(s))z(s) ds 



dy Jan dy 

+ / ^-( s >2/( s )> M ( s )) w ds ( 36 ) 
Jan ou 



'an 

where z is the solution of (fT9|) and (|20|) . 

Let w be the solution of (IMj) and (T5oT) . Subtracting w times ( TT§j) from z times 
and applying the Green's second identity yields 



f -J-(x,y)z dx = ( (wAz — zAw)dx = f (zd t 
Jn oy Jn Jan 



w — wd u z) ds 



I \ zd v w - w[-A dn z + ^(s,y,u)z + ^(s,y,u)v}\ ds 
Jan I 



dy- {s ^ u)z+ al' 



wAgnz ds 
an 



+ / <{ %w - ^0, y, u)w)z - ^(s, y, u)vw \ ds. 



an 

On the manifold dQ there holds a boundary version of Green's identity, see Lemma 
16.1 of [7]. Using this boundary version of Green's identity we have 

wAgnz ds = zAgnw ds 
nn Jan 
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and hence 
df 



dy 



x, y)z dx 



[A 9n w + d v w - ^(s, y, u)w]z - y, u)wv\ ds 

an L. dy du J 



(_~dy V ' ~~ ~du V ' u ' wv ) 



A substitution of this into (136|) gives 



<.n«),«> 



<)<) 



do 



-q-\s,V,u) - — {s,y,u)w 



v ds. 



(37) 



4 Main result 



A function u G Z7 is said to be a local solution, or a locally optimal control, of (@| 
if there is a number 5 > such that J(u) > J(n) holds for all u G Z7 satisfying 
|n — n| < 5, with their associated state y and the state constraint on y. Our main 
result is the first order necessary condition for a u G U to be a local solution. 



Theorem 4.1 // u G U is a local solution of then there exist a real number 
A > 0, a function y G C 2,a (Q) and a function w G C 2,a {pL) such that y satisfies 



— Ay + y = in Q, 

A 9n y + d v y = cp(-, y, u) on dVt 



(3f 



w satisfies 



Aw + w = —(-,y) + A^(-,y) 



(9 

Aanw + ^w = — (-,y,u)w- — (-,y,u) - A— on <9fi (39) 
<7y ay ay 



and 



/or all u E U. 



dg ( - -n dip, _ 1 



(u — u) ds > 



(40) 



Proof. Let y = G(n) be the solution of (1381) corresponding to n and p(n) = 
F(G(u)) where F is the constraint functional in our optimal control problem (j3J). 
By the theorem of Lagrange multiplier there is a A > such that 

(J'(u) +Xp'{u),u-u) > (41) 
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for all u E U . In order to obtain the necessary conditions stated in the theorem the 
only thing remaining is to calculate the derivative p'(u) of the constrain functional 
F. As in the proof of Theorem 13.21 we have 



(p\u),v) 



s, y(s))G'(u)v ds 



^(x,y(x))G'(u)v dx + [ 
oy Jan oy 

— (x, y(x))z(x) dx + / — (s,y(s))z(s)ds 



dy 



dy 



where z is the solution of (fT9|) and (|20|) corresponding to y and u. 
Let Wi be the solution of 

-Awi + Wi = — -,!/ m 12, 

oy 

Advwi + d v wi = ^-(■,y,u)w 1 - ^-(-,y,u) on dQ. 
oy oy 

By the same argument as in the proof of Theorem 13.21 we obtain 



(p'{u),v) 



dn L 



--^{s,y,u)wi 



v ds. 



(42) 



(43) 



Finally we put w = w + Xwi where w is the solution of (13"4"|) and (1331) corresponding 
to y and -u. Then u) is the solution of fl3"U|) and 



(J'(w) + Xp'(u),u- u) 



<1V 



dg 
du 



dip 
du 



(s,y,u)w 



(u — u) ds. 



5 Remarks 

A second order necessary optimality condition for the problem (T4|) can be easily 
established using exactly the same arguments as in Section [3] and HI To keep this ar- 
ticle short we leave the derivation and formulation for such a result to the interested 
readers. 

As mentioned in the Introduction the initial-boundary value problem of a parabolic 
equation with a parabolic Venttsel boundary condition arises in the engineering 
problem of heat conduction. When one considers the optimal control problem (J4]) 
with such a state equation, similar optimality conditions are expected. The study 
on such a problem is currently undergoing and the result will be published in the 
near future. Here we just point out the formulation of the state equation under the 
consideration: 

% - Ay + y = in Q, 

dt 

dy 

— - A on y- d v y = (f(-,y,u) on dtt. (44) 
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The theoretical frame work in dealing with such a problem has been well established. 
This includes existence and uniqueness of the solution in a proper function space, 
as well as the a priori estimates. For details, see [2] and [10J and the references 
therein. 

Finally we would like to point out that there is no difficulty to extend the results 
in this paper to the case when the state equation is a general second order elliptic 
equation with a general Venttsel boundary condition: 

a'Wijy + b l D { y + cy = in Q, 

a^dijy + d u y = tp(-, y, u) on dVt (45) 

where a % \ b\ c, are all constants satisfying the conditions: 

(i) {a* J } is a positive definite symmetric matrix with the smallest eigenvalue a > 

and c < 0; 

(ii) {a 1 - 7 } is a positive definite symmetric matrix with the smallest eigenvalue 7 > 0. 

In this general case, when proving a similar result to Theorem 13.21 due to the lack 
of Green's second identity, the result will take a more complicated form. 
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